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1.  Introduction 


In  a  recent  paper  a  new  method  was  described  to  evaluate  the  short-exposure  Modulation 
Transfer  Function  (MTF)  (Tofsted,  201 1).  Part  of  that  reanalysis  evaluated  diffraction-related 
effects  on  the  phase  structure  function  which  impacted  both  the  angle-of-arrival  variance  and  a 
phase-tilt  correlation  term.  However,  that  analysis  considered  only  Kolmogorov  turbulence  (no 
outer  scale  of  turbulence  influence).  The  outer  scale  is  well  known  to  affect  angle-of-arrival  (e.g., 
Tofsted,  1992),  but  that  and  other  analyses  did  not  consider  aperture  effects  leading  to  diffraction 
influences.  To  study  the  outer  scale  influence,  the  low  frequency  portion  of  the  refractive  index 
spectrum  must  first  be  modeled,  then  its  impacts  on  the  angle-of-arrival  must  be  evaluated 
numerically.  This  low  frequency  influence  is  parameterized  relative  to  the  outer  scale  of 
turbulence,  Lo. 

For  completeness,  we  begin  by  considering  prior  angle-of-arrival  variance  formulations,  as 
reported  by  Beland  (1993)  and  Dror  and  Kopeika  (1995).  Beland  (1993)  provided  a  formula  for 
plane  wave  angle-of-arrival  variance: 

(a^)  =  2.91  r  Cl{z)  d;^,  (1) 

Jo 

where  D  is  the  system  entrance  pupil  diameter  [m],  L  is  the  path  length  [m],  2:  is  the  path 
variable,  and  is  the  refractive  index  structure  parameter  [m“^/^]. 

Dror  and  Kopeika  (1995)  suggested  the  spherical  wave  single-axis  angle-of-arrival  variance 
formula, 

{al)  =  2.92  C^^iz)  (^-j  d^,  Z  <  1  «  Q,  (2) 

where  diffraction  and  outer  scale  related  dimensionless  scales  are  introduced:  Q  =  D/{\ L)^/^ 
and  Z  =  D/Lo,  where  A  is  the  propagation  wavelength.  (Note,  Dror  and  Kopeika  claimed 
equation  2  was  a  single-axis  result,  but  comparison  with  equation  1  indicates  a  discrepancy. 

Either  the  Dror  and  Kopeika  result  is  too  large  by  a  factor  of  two,  or  the  Beland  result  is  too  small 
by  a  factor  of  two.) 

The  current  study  extends  the  above  results  to  include  outer  scale  and  diffraction  effects 
(dependence  on  Q  and  Z  variables)  for  the  spherical  wave  case.  For  tabulation  purposes  the 
diffraction  and  outer  scale  effects  are  studied  under  conditions  of  homogeneous  turbulence. 
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Motivating  this  study  is  the  application  where  measured  angle-of-arrival  variance  is  used  to 
predict  turbulence  strength,  That  is,  experiments  specifically  designed  to  study  the  impacts 
of  turbulence  typically  measure  turbulence  strength  directly  using  one  or  more  scintillometers.  In 
other  circumstances  collected  imagery  may  exhibit  turbulence  impacts  but  lack  corresponding 
turbulence  measurements.  Further,  some  imaging  configurations  may  preclude  the  use  of 
scintillometers,  for  example  if  objects  of  interest  are  airborne.  Hence,  if  imagery  alone  can  be 
analyzed  to  estimate  this  extends  the  potential  analysis  of  turbulence  effects  to  greater 
numbers  of  data  sets. 

For  this  case,  it  is  important  to  know  how  outer  scale  and  diffraction  effects  are  impacting  the 
estimates  of  turbulence  via  angle-of-arrival.  As  will  be  shown,  while  the  error  is  not  great, 
including  outer  scale  effects  adds  an  additional  source  of  error  to  the  calculations.  The  remainder 
of  the  paper  is  devoted  to  resolving  how  to  evaluate  the  angle-of-arrival  variance  given  knowledge 
of  the  turbulence  strength,  outer  scale,  and  optical  problem  geometry.  We  first  overview  the 
calculation  in  section  2,  but  leave  the  final  resolution  of  the  analysis  to  section  4.  Between  these 
two,  in  section  3  we  introduce  the  revised  von-Karman  (RvK)  spectrum  first  described  in  Tofsted 
et  al.  (2007),  including  the  definition  for  the  outer  scale  given  in  terms  of  height  above  the 
ground  and  the  Obukhov  length  (e.g.,  Obukhov,  1946;  Paulson,  1970). 


2.  General  Angle-of-Arrival  Variance  Calculation  Procedure 


In  this  section  we  overview  the  angle-of-arrival  variance  calculation  procedure.  Let  be  the 
angle-of-arrival  variance,  assumed  isotropic  such  that  x  and  y  axis  components  are  equal.  To 
model  this  variance,  the  basic  tilt  calculation  included  in  Fried’s  (1966)  short-exposure  imaging 
paper  is  adapted.  That  approach  considered  the  phase  produced  by  a  point  source  and  propagated 
a  distance  L  to  the  entrance  pupil  of  a  circular  aperture  system.  There,  Fried’s  main  integral  form 
was  integrated  in  the  aperture.  Fried  considered  the  propagated  phase  could  be  expressed  as  a 
function  0('u),  of  the  two-dimensional  vector  in  the  transverse  plane  of  the  receiver  entrance 
pupil,  V. 

Fried  considered  that  the  tilt  effect  could  be  parameterized  by  a  vector,  a,  whose  dimensions  are 
that  of  inverse  length  (m“^),  to  represent  the  mean  tilt  of  the  incident  phase  across  the  aperture. 
This  vector  permits  a  new  tilt-corrected  phase  function  to  be  defined,  as  ijj{v)  =  0(^7)  —  a  ■  v,  in 
the  system  entrance  pupil. 
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The  variance  of  a  could  be  expressed  in  integral  form  in  the  system  aperture  by, 

=  7=;2  /  [^ciu)  -  Tl{u)]  V^{Du),  (3) 

L>^  Jo 

where  functions  JFc{u)  and  JFl{u)  are  given  by, 

^c{u)  =  (2/7r)  (^cos“^(ci;)  —  a;  Vl  —  , 

J^l{u)  =  (2/7r)  (^3  cos“^(ci;)  —  a;  7  — Vl  —  ,  (4) 

and  where  'D^{p)  is  the  phase  structure  function  (e.g.,  Goodman,  1985)  for  a  transverse  separation 
distance  p  in  the  system  aperture. 

The  vector  a  physically  represents  the  rate  of  phase  change  per  unit  distance  in  the  system 
aperture.  This  vector  is  transformed  into  an  angular  dimension  by  dividing  the  variance  of  a  by 
the  wavenumber  {k  =  2ti /X).  The  a  variance  can  thus  be  transformed  into  an  angular  variance 
via: 

{ct  ■  cx)  =  {a  ■  a) /k"^ .  (5) 

This  transformation  is  based  on  an  argument  using  similar  triangles. 

The  key  to  extending  these  results  beyond  Kolmogorov  turbulence  is  to  include  outer  scale 
influences  in  the  phase  structure  function  calculation.  For  spherical  wave  propagation  this 
function  is  given  as, 

VJp)  =  k"^  L  [  dc  [  KdK^niK-)  [1  —  Jo{kcp)]  cos^  - ^  .  (6) 

Jo  Jo  [  47r 

where  c  is  a  normalized  path  variable  (c  =  z/L),  k  [m“^]  is  the  spatial  frequency  in  which  the 
turbulence  spectrum  $„(«;)  is  expressed,  and  Jo{x)  is  the  Bessel  function  of  the  first  kind  order 
zero. 

Notice  that  in  equation  3  utilizes  the  normalized  distance  argument  u  =  p/D.  In  Tofsted 
(201 1),  it  was  shown  that  when  invoking  constant  Kolmogorov  turbulence  the  phase  structure 
function  could  be  expressed  as  a  product  of  the  wave  structure  function,  'P(p),  with  a  function, 
denoted  a{R),  accounting  for  the  diffraction  influence  arising  from  the  argument  of  the 
cosine- squared  term  in  the  phase  structure  function’s  inner  integral.  Comparison  with  the  wave 
structure  function’s  integral  definition  is  useful: 

p1  pod 

V{p)  =  k"^  L  dc  KdK^nin)  [1  —  Jo{kc p)]  .  (7) 

Jo  Jo 
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The  phase  structure  function’s  integral  form  approaches  that  of  the  wave  structure  function 
whenever  the  argument  of  the  cosine-squared  term  remains  small  throughout  integration.  The 
argument  of  a{R)  is  the  normalized  length  scale  R  =  p/P  where  P  =  y/XL,  the  Fresnel  zone, 
such  that  R  =  u  D/ P  =  uQ. 

This  construction,  V^{uD)  =  'D{uD)  a{uQ),  is  possible  because  the  Kolmogorov  refractive 
index  spectrum, 

^n{n)  =  ^K{n)=(3ClK-^^/\  (8) 

contains  only  the  single  parameter,  and  where  /9  is  an  integration  constant. 


13  =  (5/36)  [22/3r(5/6)]  / 


'7r3/2r(2/3) 


^  0.033. 


(9) 


connecting  the  Kolmogorov  spectrum  to  the  unrestricted  refractive  index  structure  function. 


Vn{r)  =  Vxir)  =  Clr‘^1^. 


(10) 


In  the  case  of  homogeneous  turbulence,  C\  factors  out  of  both  structure  function  integrals.  In  the 
case  of  the  phase  structure  function,  the  remaining  double  integral  becomes  a  function  of  only  R. 

This  simplified  result  is  considerably  complicated  by  the  introduction  of  an  outer-scale 
dependence  in  the  refractive  index  spectrum,  as  discussed  in  the  following  section. 


3.  Revised  von  Karman  Spectrum 


The  influence  of  the  outer  scale  of  turbulence  on  angle-of-arrival  measurements  is  to  introduce  a 
low  frequency  cutoff  in  the  spectrum.  The  form  used  to  model  the  outer  scale  effect  was  first 
introduced  in  2007  (Tofsted,  2007).  The  revised  von-Karman  (RvK)  spectrum  is  based  on 
von-Karman’s  (1948)  recommendation  that  the  low  frequency  dependence  in  the  spectrum  should 
exhibit  behavior.  In  this  section  we  describe  this  spectrum  and  define  the  outer  scale  behavior 
in  the  atmospheric  surface  layer.  In  particular,  we  define  the  outer  scale  as  a  function  of  height 
above  the  surface  (optical  line-of-sight  height)  as  well  as  the  Obukhov  length  Lqb  that 
characterizes  vertical  structure  within  the  surface  layer.  This  analysis  was  originally  discussed  in 
Tofsted  (2000)  and  is  based  on  analysis  by  Kaimal  et  al.  (1972)  of  data  collected  during  the  1968 
Kansas  experiment.  Due  to  space  considerations,  the  interested  reader  is  directed  to  the 
referenced  documents,  since  only  the  equations  describing  the  spectrum  and  outer  scale  are 
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provided  here. 


The  RvK  spectrum  can  be  expressed  as  the  product  of  a  Kolmogorov  spectrum,  from 

equation  8,  and  a  high-pass  filter  component, 

Lo)  =  ^k{k)^{kLo),  (11) 


where  the  outer  scale  modulating  term  is  expressed  using  a  weighted  sum, 

i=3 

^{x)  ='^WiUi{aix),  (12) 

i=l 


with  weighting  factors,  Wi  =  (0.25,  0.50,  0.25)  and  frequency  coefficients 
tti  =  (0.80,  2.00,  5.00),  respectively,  in  the  basic  high-pass  function, 

U,{X)  = 

(1  +  X2)^^/® 


(13) 


Overall,  the  RvK  spectrum  is  consistent  with  the  spectrum  of  Tofsted  (2007),  but  the  above 
derivation  clearly  shows  the  filtering  effect  of  Lo).  The  Ui  function  also  guarantees  that  the 
spectral  dependence  at  low  frequency  is  consistent  with  von-Karman’s  low  frequency 
prescription,  and  with  the  required  inertial  subrange  high  frequency  dependence. 


The  rationale  for  employing  three  shifted  copies  ofUi  is  based  on  comparison  between  the  form 
of  Ui  and  the  plots  of  Kaimal  et  al.  (1972).  The  Ui  function  is  narrower  than  reality,  requiring 
multiple  shifted  copies. 


At  this  point  perhaps  it  would  be  valuable  to  make  some  comment  with  regard  to  the  traditional 
form  taken  by  the  so-called  von-Karman  (vK)  spectrum.  First,  how  this  spectrum  came  about 
and  who  originally  proposed  this  form  is  somewhat  of  a  mystery,  though  it  does  appear  in  early 
writings  of  Tatarskii  (1971).  It  is  commonly  expressed  as. 


Lo)  —  (3 


1 

(L-2  +  fi;2)ll/6 


(14) 


Note  that  this  spectrum  reaches  a  maximum  at  zero  frequency,  while  von-Karman  specifically 
indicated  a  low  frequency  behavior  of  Also,  recall  that  the  refractive  index  spectrum  is 
actually  the  power  spectrum  for  a  space  containing  refractive  index  perturbations.  That  is,  one 
subtracts  from  the  spatially  varying  refractive  index  function  its  mean  value.  Therefore  the  space 
of  refractive  index  perturbations  has  a  zero  mean.  But  the  zero  frequency  component  of  the 


5 


power  spectrum  is  just  the  squared  mean  value,  which  is  zero.  Therefore  the  traditionally  used 
vK  spectrum  is  technically  incorrect  while  the  RvK  spectrum  properly  characterizes  the  behavior 
at  the  origin  and  reflects  an  energy  containing  range  (e.g.,  Hinze,  1987).  This  is  not  to  say  that 
this  spectrum  is  without  faults.  It  assumes  isotropy,  which  is  inaccurate  for  larger  scale 
perturbations,  but  that  is  a  topic  for  a  separate  analysis.  Perhaps  a  more  accurate  spectrum  could 
be  developed  in  the  future  that  includes  a  flattening  of  the  spectrum  at  lower  frequencies, 
consistent  with  more  recent  analyses  (e.g.,  Tofsted  et  ah,  2009). 

Next,  we  consider  parameterizing  the  outer  scale.  KaimaTs  results  also  included  data  that  could 
be  analyzed  to  form  a  relationship  between  the  outer  scale  and  height  above  the  surface,  H,  and 
the  Obukhov  length,  Lqb-  In  this  process  it  also  became  obvious  that  standard  spectral  forms 
containing  an  outer  scale  simply  included  Lo  as  a  parameter,  but  without  attaching  that  parameter 
to  any  physically  meaningful  property  of  the  atmosphere.  As  such,  the  best  course  appeared  to 
argue  that  a  specific  definition  of  outer  scale  be  stated:  Let  Ko  =  l/Lo,  such  that  the  outer  scale 
may  be  defined  as, 

Lo) I =  1/2.  (15) 

The  outer  scale  is  thus  just  the  3-dB  level  of  the  Revised  von-Karman  spectrum  relative  to  the 
Kolmogorov  spectrum.  As  a  matter  of  interest,  one  may  also  define  an  integral  scale,  Ao, 
representing  the  peak  of  the  overall  spectrum,  which  we  expect  to  occur  roughly  in  the  center  of 
the  “Energy  Containing  Range”  (Hinze,  1987).  This  scale  occurs  approximately  at  Ao  =  2  Lq. 

Based  on  this  outer  scale  definition,  the  Kaimal  et  al.  plot  results  could  be  analyzed  (Tofsted, 
2000),  and  the  outer  scale  modeled  as. 


Lo  = 


(77/(4.34  7+^  +  0.013),  ^  >  0, 

1  77/(0.34  7^  +  0.013),  e<0. 


(16) 


where  ^  =  Hj Cob  is  a  dimensionless  height  ratio.  Note  that  for  large  ^  argument,  the  outer  scale 
increases  roughly  as  a  function  of  (77  \  Cob\)^^‘^-  In  general  Cos  <  0  during  daylight  hours  when 
the  sensible  heat  flux  is  positive  (ground  is  heating  the  air),  and  Cqb  >  0  at  night  when  this 
situation  is  reversed.  Thus  Lo  is  generally  larger  during  the  daytime  than  at  night,  even  at  the 
same  height,  a  result  that  generally  conflicts  with  the  commonly  used  approximation  Lo  ^  kH, 
where  k  ~  0.4  is  von-Karman’s  constant. 


To  determine  the  impact  of  this  spectrum  on  angle-of-arrival,  we  examine  the  influence  of  the 
spectrum  on  wave  and  phase  structure  functions,  and  the  resulting  influence  on  the 
angle-of-arrival  variance  based  on  the  equations  in  section  2.  In  order  to  compare  these  new 
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results  with  the  prior  Kolmogorov-based  results,  the  new  expressions  will  be  derived  in  forms 
consistent  with  Fried’s  original  analysis. 

4.  Outer  Scale  and  Diffraction  Influenced  Angle-of-Arrival 

Using  the  Fried  angle-of-arrival  formulation  described  in  section  2,  we  first  describe  the  results  of 
that  analysis  for  Kolmogorov  turbulence,  then  consider  the  influence  of  the  outer  scale  along  with 
diffraction  effects.  However,  from  the  point  of  view  of  developing  the  results,  we  need  to  work 
backwards,  starting  with  the  wave  structure  function,  then  dealing  with  the  phase  structure 
function,  which  contains  a  diffraction  related  component,  and  finally  the  Kolmogorov  based 
angle-of-arrival  variance.  We  can  then  work  this  same  process  for  the  Revised  von-Karman 
spectrum. 

Let  us  begin  by  repeating  the  equation  for  the  wave  structure  function: 

p1  roo 

V{p)=8n'^k'^L  dc  KdK^n{n)  [^  —  Jo{ncp)]  .  (17) 

Jo  Jo 

To  simplify  the  integrals,  we  let  U  =  ncp  =  kcDu  replace  the  Bessel  function  argument.  We 
then  transform  the  argument  of  the  spectral  inner  integral  to  one  extending  over  V : 

/*!  foo 

V(Du)  =  [i  _  Jo(\/)]  dU.  (18) 

Jo  Jo 

For  constant  turbulence  conditions  the  inner  and  outer  integrals  are  separable,  where, 

rl  Q  roo 

/  c+5/3dc=-,  /  -  Jo(U)]  dU  =  K  ^  1.1183344.  (19) 

Jo  8  Jo 

These  results  allow  us  to  define  a  nominal  coherence  length, 

p;^/^  =  4x^71^  (J^k^LCl  (20) 

O 

such  that  the  wave  structure  function  for  constant  Kolmogorov  turbulence  may  be  expressed, 
using  S  =  D/po, 

V{Du)  =  2  ^+5/3  ^  2  ^5/3  ^5/3_  (21) 

Using  this  form,  we  next  consider  the  effect  of  diffraction  on  the  phase  structure  function.  This 
result  was  previously  discussed  in  Tofsted  (201 1).  As  before,  we  repeat  the  phase  structure 
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function  equation: 


p1  poo 

VJp)  =  STi‘^k‘^L  dc  ndn^niK-)  [1  — </o(Kcp)]  cos^ 
Jo  Jo 


XLc{l  —  c) 


47r 


(22) 


Again  introducing  the  same  dimensionless  variable  V  =  ncp  =  Duck,  and  again  utilizing  po,  the 
double  integral  form  can  be  expressed  as, 


V^{Du) 


2 


Dc^/^dc  1/-8/3 

AT  ;  Jo  TW  D 


[1  -  Jo{V)]  cos^ 


1/2  c(l  -  c) 
{QuY  dvrc^ 


(23) 


Clearly  the  double  integral  is  only  a  function  of  Qu,  reflecting  the  influence  of  diffraction. 
Following  a  similar  methodology  as  in  previous  papers  on  the  topic  of  the  phase  structure 
function,  we  now  write  this  function  as. 


V^{Du)  =  2  (^m)^/3  a{Qu)  =  V{Du)  a{Qu), 


(24) 


where. 


a{Qu) 


C  dc 

0  TW 


1/-8/3 


[1  -  Jo{V)]  cos2 


1/2  C(1  -  C) 

{Qu^  dvrc^ 


(25) 


The  form  taken  by  this  function  is  plotted  in  figure  1  for  ten  decades  of  R  =  Qu  variation. 


Figure  1 .  Phase  structure  function  diffraction  dependence  factor 
a{R),  with  R  =  Qu,  for  constant  Kolmogorov 
turbulence. 
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Introducing  the  equation  24  representation  of  the  phase  strueture  funetion  into  the  equation  3  tilt 
parameter  varianee  equation,  and  then  introdueing  this  result  into  equation  5,  we  produee  an 
equation  for  the  angle-of-arrival  varianee  of 

1 28  /*! 

lo?)  =  p  ^2  du  [J^c{u)  -  J^l{u)\  u^/^a{Qu).  (26) 

The  integral  over  the  system  aperture  is  elearly  a  funetion  of  Q  alone,  whieh  may  be  further 
normalized  relative  to  the  limiting  response  at  Q  — C)0,  where  a  =  1.  In  this  limit,  the 
eomponent  integrals  evaluate  to, 

duTc{u)  ^  0.03749;  duTiiu)  ^  0.00489.  (27) 

Jo  Jo 

sueh  that  128  x  (0.03749  —  0.00489)  =  4.1728,  and  the  angle-of-arrival  varianee  may  be 
expressed  as, 

(a^)  =  4.1728  G{Q)/{k^  D^).  (28) 

The  funetion  G{Q),  thus  dehned,  is  plotted  in  hgure  2.  This  funetion  assumes  near-limiting 
behavior  G{Q)  1  for  Q  >  3. 


Figure  2.  Angle-of-arrival  varianee  diffraction  dependence 

function  G{Q),  with  Q  =  i7/\/XL  =  D/P,  for  constant 
Kolmogorov  turbulence. 


Comparison  of  this  result  with  the  Dror  and  Kopeika  (1995)  and  Beland  (1993)  equations  requires 
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expanding  =  {D/pof’^^  such  that  the  limit  at  large  Q  can  be  expressed  as, 

=  4.1728  (3/2)  7r2/3KLC'2/Di/3  =  2.2798 LC'2/D1/3.  (29) 

For  constant  Dror  and  Kopeika’s  equation  leads  to  a  scaling  factor  of  2  x  2.92  x  (3/8)  =  2.19 
corresponding  to  the  full  angle-of-arrival  variance.  Comparing  the  lead  constants  of  Beland 
(equation  1)  and  that  used  by  Dror  and  Kopeika  (equation  2),  apparently  the  Beland  result 
corresponds  to  a  single-axis  angle-of-arrival,  rather  than  the  full  angle-of-arrival  variance. 

In  a  more  recent  publication  (Tofsted,  2012)  the  question  of  path  varying  turbulence  strength  was 
addressed,  revealing  that  this  variation  can  be  accounted  for  through  rescaling  of  the  S  and  Q 
parameters.  Thus  the  above  forms  are  identical  for  variable  path  height  geometries.  However, 
the  introduction  of  path  varying  outer  scale  might  introduce  unusual  changes  in  the  wave  and 
phase  structure  functions.  Therefore,  in  this  initial  study,  only  path  constant  outer  scale 
influences  will  be  considered. 

Given  that  caveat,  we  are  now  prepared  to  extend  the  angle-of-arrival  analysis  to  the  case  of  the 
RvK  turbulence  spectrum.  Following  the  pattern  established  when  studying  angle-of-arrival  for 
Kolmogorov  turbulence,  we  first  direct  our  attention  on  the  wave  structure  function.  Introducing 
the  outer  scale  influenced  spectrum  in  equation  17,  the  argument  to  the  Bessel  function  can  again 
be  replaced  by  the  variable  V,  but  the  resulting  double  integral  no  longer  factors  as  two 
independent  integrals.  Instead,  a  dependence  on  the  dimensionless  variable  Zuc  appears: 

V(Du)  =8n^  (3eLCl{DuY^/^  c+^/^dc  dH.  (30) 

Jo  Jo  \ZucJ  F®/'* 

We  may,  nevertheless,  divide  through  the  double  integral  by  the  factor  (3/8)  K  such  that  we  may 
again  express  the  function  in  terms  of  the  nominal  coherence  length  po  from  equation  20.  The 
complete  wave  structure  function  may  thus  be  expressed  using  a  form  similar  to  equation  21, 

V{Du)  =2{Su)^/^  H{Zu),  (31) 


where  the  new  outer-scale  dependent  function. 


H{Zu)  = 


1  c+5/3 


dc 


Q 


V  \  [1- Jo(V^)] 


ZucJ  KH8/3 


dH, 


(32) 


is  introduced.  Now  note  that  the  D  factor  approaches  unity  at  large  argument.  Thus,  for  small  Z, 
the  integral  approaches  the  limiting  Kolmogorov  behavior.  That  is,  H (Zu)  must  have  a  limit  of 
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unity  at  small  Zu.  Conversely,  for  large  Z,  the  outer  scale  influence  will  dominate.  In  this  limit 
it  is  advantageous  to  introduce  yet  another  dimensionless  constant,  M  =  Lo/po,  along  with  the 
function  Hh{x)  =  H{x),  such  that, 

V{Du)=  2  Hh{Zu).  (33) 

Hh{Zu)  approaches  a  constant  limit  as  Zu  oo,  V{oo)  ^  8.40  Hh{oo)  =  4.20.  The 
function  Hh{Zu)  is  plotted  in  figure  3,  along  with  the  small  scale  asymptotic  dependence 

{Zuf!^. 


Figure  3.  Outer  scale  cutoff  of  the  wave  structure  function  (solid 
line)  compared  to  power  law  for  Kolmogorov 

wave  structure  function. 


The  effect  of  the  outer  scale  limitation  on  the  refractive  index  spectrum  is  such  that  the  wave 
structure  function  exhibits  an  upper  limit.  Equivalently,  the  function  H{Zu)  must  be 
proportional  to  {Zu)~^/^  for  large  Zu  such  that  'D{p)  approaches  a  finite  upper  bound. 

Using  the  modified  equation  33  wave  structure  function,  we  next  explore  the  extension  of  the 
outer  scale  impacts  to  the  phase  structure  function.  To  express  this  structure  function  for  the  RvK 
model,  we  again  use  the  V  variable,  and  pull  the  cosine-squared  term  from  equation  25,  to 
produce. 


V^{Du) 


2  (^m)'^/3 


^1^5/3  dc  ^  ^ 

Jo  {3/S)  Jo  K  1/8/3  ^\zuc) 


1/2  (1  -  c) 

(Qm)2  Attc 


(34) 
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Following  the  model  of  the  analysis  for  Kolmogorov  turbulenee,  the  wave  strueture  funetion 
effeets  are  extraeted,  leaving  a  funetion  of  both  Qu  and  Zu,  refleeting  the  influenees  of  both 
diffraetion  and  outer  seale.  The  phase  strueture  funetion  may  thus  be  expressed  as, 

V^{Du)  =  V{Du)  A{Qu,  Zu)  =  2  (Su)^^^  H{Zu)  A{Qu,  Zu),  (35) 


where, 


A{Qu,  Zu)  =  [ 
Jo 


1  dc 


[1- Jo(f")]  n[v/{zuc)] 
K  1/8/3  H{Zu) 


1/2  (1  -  C) 

(Qu)'^  Att  c 


(36) 


The  form  taken  by  this  two-dimensional  funetion  is  plotted  for  a  series  of  values  of  Qu  as 
funetions  of  Zu  in  figure  4.  Zu  varies  over  the  range  from  10“®  to  exhibiting 
near-asymptotie  behaviors  at  both  extremes.  The  Qu  values  ranged  from  Qu  =  2“®,  for  the 
lowest  curve,  to  2+®  for  the  highest  curve,  in  increments  of  factors  of  2  between  consecutive 
curves. 


Figure  4.  Outer  scale  impact  on  transition  function  of  phase 

structure  function.  13  plots  of  different  Qu  values  from 
(bottom-most  curve)  to  (upper-most  curve) 
differing  by  a  factor  of  2  between  each. 


To  compare  the  behavior  of  A{Qu,  Zu)  with  that  of  a{Qu),  we  consider  the  limit  Zu  <C  1.  To 
compare  these  cases  we  plot  the  A{Qu,  10“®)  against  a{Qu)  in  figure  5.  As  expected,  the 
asymptotic  behavior  of  A{Qu,  Zu)  ranges  between  0.5  and  1.0  at  Zu  -C  1,  corresponding  to  the 
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asymptotic  behaviors  of  Qu  -C  1  and  Qu  S>  1,  respectively,  seen  in  figure  1.  The  range  of  Qu 
selected  is  characteristic  of  most  imaging  scenarios  except  for  extremely  large  or  small  optics. 


logio  (R) 


Figure  5.  Comparison  of  a{R)  from  the  Kolmogorov  turbulence 
model,  versus  A{R,  Zu)  from  the  Revised  von  Karman 
turbulence  model,  with  R  =  Qu. 


For  the  remaining  behavior  of  A{Qu,  Zu),  obviously  the  function  exhibits  a  complicated  behavior 
that  defies  easy  characterization.  So,  rather  than  attempting  to  characterize  this  intermediate 
result,  let  us  instead  continue  on  to  determine  the  angle-of-arrival  variance  for  RvK  turbulence. 

Introducing  the  tabulated  A{Qu,  Zu)  results  into  a  numerical  integration  routine,  we  begin  with 
equation  26  as  modified  by  the  RvK  phase  structure  function: 

1  no  c5/3  ^1 

{a^)  =  p  ^2  Jq  [^c{u)  -  J^l{u)]  u^^^  H{Zu)  A{Qu,  Zu).  (37) 

Now,  the  results  of  this  integration  will  approach  those  of  the  equation  27  limiting  behaviors 
when  Zu  1  -C  Qu.  Therefore,  for  comparison  purposes  we  may  multiply  and  divide  by  the 
factor  0.03260  =  0.03749  -  0.00489.  This  produces  the  same  lead  coefficient  as  in  equation  29 
such  that  we  can  write, 

(a")  =  2.2798  L  Cl/ D^/^Gg{Q,Z),  (38) 

where, 

Gc(Q,  Z)  =  /'  dt,  H(Zu)  A(Qu.  Zu).  (39) 

JO  0.03260 
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Like  the  G{Q)  funetion,  the  Gc{Q,  Z)  funetion  will  approaeh  unity  when  Z  -C  1  -C  Q,  but  will 
then  deeay  when  this  eondition  does  not  hold.  This  funetion  is  plotted  for  a  range  of  Q  values  as 
a  funetion  of  Z  in  figure  6. 


Figure  6.  Weighting  function  Gg{Q,  Z)  for  the  range 

10“®  <Z  <  10+^  and  2“®  <  Q  <  2+®  (lowest  to 
highest  curves). 


Figure  7  then  shows  the  similarity  in  the  limiting  behavior  of  GciQ,  Z)  at  low  values  of  Z 
eompared  to  the  behavior  of  G{Q)  from  figure  2.  From  this  eomparison  we  observe  an  apparent 
eorrelation  between  the  effeets  of  small  Z  and  small  Q  in  the  Gg{Q-,  Z)  model  sueh  that  at  low  Q 
Gg  traeks  inereasingly  below  G{Q)  as  Q  drops. 


5.  Discussion 


We  thus  have  developed  a  pieture  of  the  behavior  of  the  angle-of-arrival  varianee  as  the  outer 
seale  and  diffraetion  effeets  vary.  However,  as  we  also  have  a  means  of  eonneeting  the 
atmospherie  propagation  eharaeteristies  to  typieal  meteorologieal  eharaeteristies,  let  us  eonsider 
what  we  might  eneounter  as  typieal  effeets.  To  begin,  let  us  seleet  a  typieal  windspeed  of 
3.5  m/s.  This  is  eharaeteristie  of  a  daytime  desert  atmosphere  (e.f.,  figure  14  of  Tofsted  et 
al.  (2009),  where  we  eonsidered  mean  hourly  windspeeds  eolleeted  over  a  60  day  period  in  2008). 
Let  us  also  eonsider  a  path  height  of  =  4  m  AGL  (above  ground  level).  Typieal  daytime 
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Figure  7.  Comparison  of  G{Q)  function  plotted  in  figure  2  against 
the  Gg{Q,  Z)  from  figure  6  for  small  Z. 


temperatures  during  our  measurement  period  averaged  around  25  °C  or  roughly  300  K.  The 
White  Sands  area  stands  at  approximately  1400  m  ASL  (above  sea  level)  with  typical  atmospheric 
pressure  of  P  =  880  mbar  (millibars),  and  characteristic  sensible  heat  flux  values  were  in  the 
range  200-600  W/m^.  We  write  the  equation  for  the  Obukhov  length  as  (Tofsted,  1993), 


,  pCpTul 

J^Ob  — - j - TO — 

kg  Us 

where  p  is  the  dry  air  density, 

p=  0.34838  P/T, 


(40) 


(41) 


with  pressure  P  in  millibars,  and  temperature  T  in  Kelvin,  such  that  p  is  expressed  as  kg/m^;  Cp 
is  the  heat  capacity  of  air  at  constant  pressure,  1004.84  J/(kg  K);  k  is  von-Karman’s  constant, 
k  ~  0.4;  g  is  the  gravitational  attraction,  g  =  9.81  m/s^;  m*  is  the  friction  velocity  in  m/s;  and  1-Ls 
is  the  sensible  heat  flux.  While  the  true  friction  velocity  exhibits  an  effect  due  to  diabatic 
influences,  for  the  sake  of  simplicity  here  we  use  only  the  logarithmic  wind  profile,  such  that  ijjm, 
the  momentum  diabatic  function,  is  zero,  and  we  posit. 


ku  ku 

ln(P/Po)  -  ~  ln(P/Po)’ 


(42) 


where  Ho  is  the  surface  roughness  length  (Hansen,  1980),  which  we  determined  through 
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measurements  to  be  approximate  4.5  cm,  which  approximately  matches  Hansen’s  recommended 
relative  length  compared  to  typical  roughness  element  heights  of  1  m  (mesquite  bushes,  grass,  and 
yuccas  plants).  For  a  2-m-wind  measurement  height,  and  wind  speed  u  =  3.5  m/s,  we  thus 
estimate  a  friction  velocity  of  m*  0.37  m/s. 

The  atmospheric  density  is  approximately  p  =  1.022  kg/m^.  Hence,  pCp  =  1026.9  J/(m^  K), 
such  that  Q  =  Hs/ipCp)  will  range  from  0.195  to  0.584  K  m/s.  Therefore,  we  can  express  the 
range  of  the  Obukhov  length  as. 

Tv? 

Lob  =  --j — ^  = -20  m...  -  7  m.  (43) 

kg  Q 

For  this  range  of  Obukhov  length  we  obtain  from  equation  16  a  range  of  outer  scale  from  14  m  to 
24  m  for  a  path  height  of  4  m,  where  the  outer  scale  is  shorter  under  stronger  turbulence 
conditions. 

Consider  then  a  range  of  optics  aperture  sizes  ranging  from  2.5  cm  for  binocular  optics  to  20  cm 
for  a  telescope  optic.  For  this  range  of  the  parameter  D,  let  us  also  consider  path  length  ranging 
from  1  km  to  4  km  and  wavelength  equal  to  0.55  pm  (visible  radiation),  or  a  range  of 
P  =  (A  L)^/^  of  2.35  cm  to  4.7  cm. 

For  D  =  2.5  cm,  Z  =  D/Lo  ranges  from  0.0010  to  0.0018,  while  Q  =  D/P  varies  from  0.53  to 
1.06.  Hence,  from  figure  6  the  Go  coefficient  varies  from  0.85  to  0.90  where  the  impacts  of  Q 
and  Z  are  approximately  equal. 

For  D  =  20  cm,  Z  ranges  from  0.0083  to  0.0143,  while  Q  varies  from  4.26  to  8.51,  and  Go  varies 
from  0.88  to  0.93  where  the  impact  is  virtually  all  due  to  variations  in  Z. 

Overall,  the  scaling  constant  from  figure  6  is  between  0.80  and  0.95,  or  (in  combination  with  the 
scaling  constant  2.28  from  equation  38)  when  compared  to  the  constant  2.19  derived  from 
equation  2,  one  obtains  a  range  from  1.824  to  2.17.  That  is,  angle-of-arrival  variance  ratios  range 
from  0.833  to  0.989.  Of  course,  if  we  intend  to  use  angle-of-arrival  variance  to  compute  G'/,  we 
recognize  that  turbulence  strength  is  highly  variable,  often  exhibiting  a  third  of  an  order  of 
magnitude  variance  around  a  mean  with  correlation  times  on  the  order  of  several  seconds  or  less. 
Also,  since  most  imaging  scenarios  involve  Q  values  greater  than  unity,  and  usually  of  the  order  2 
or  greater,  we  expect  that  the  outer  scale  influence  (rather  than  diffraction)  will  be  the  dominant 
term,  and  this  impact  does  extend  down  to  produce  a  10%  effect  even  for  Z  values  on  the  order  of 
0.01,  and  not,  as  Dror  and  Kopeika  indicated,  merely  Z  <1. 
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6.  Conclusions 


In  this  study  a  new  analysis  of  the  influenee  of  both  diffraetion  and  outer  seale  of  turbulenee 
effeets  on  angle-of-arrival  varianee  has  been  performed.  This  analysis  has  shown  that  outer  seale 
influenees  ean  be  signifieant  even  for  aperture  sizes  far  less  than  the  outer  seale  size.  To  use  the 
information  developed  in  this  report,  one  ean  simply  apply  rule-of-thumb  estimates  in 
determining  the  Obukhov  length,  plus  estimates  of  the  optieal  path  height,  to  determine  the  outer 
seale.  From  this  knowledge,  and  knowledge  of  the  system  aperture  size,  one  ean  easily  determine 
the  eorreetion  neeessary  to  obtain  a  more  aeeurate  estimate  of  from  observed  edge  motion 
analysis. 

This  analysis  has  also  shown  the  relative  strengths  and  weaknesses  in  previous  analyses.  In 
partieular,  the  Dror  and  Kopeika  (1995)  equation  appears  to  be  roughly  eorreet  as  long  as 
Z  <^1  <  Q.  However,  this  eondition  differs  from  that  indieated  by  Dror  and  Kopeika  of 
Z  <  1  Q.  More  generally,  the  new  formulation  provides  a  method  of  estimation  under  a  far 
wider  set  of  measurement  eonditions  than  that  supported  by  previous  analyses.  Also,  Tofsted 
(1992)  only  eonsidered  the  outer  edge  of  the  wave  propagating  toward  the  aperture  when 
evaluating  angle-of-arrival  varianee.  The  eurrent  analysis  is  obviously  more  suitable  for  eireular 
aperture  systems  involving  an  ensemble  of  point  separations  in  the  aperture. 
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